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>■ ' We calculate the squarks one-loop corrections to the cross sections of the charged Higgs 

^ I boson pair production in photon-photon collisions in the minimal supersymmetric extension 
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0> I squark and the right-top squark and for tan/? near 1.5, the corrections can reduce the cross 

Q^l sections by more than 10% for a wide range of the charged Higgs boson mass, depending on 
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percent. Those corrections can be comparable to the 0{am^/m^) Yukawa corrections in 

/S ' previous literature. 
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I. INTRODUCTION 

In various extension of the Higgs sector of the standard model (SM), there are physical 
charged Higgs bosons. For example, the minimal supersymmetric standard model(MSSM) 
|l[] with two Higgs doublets, which give masses separately to up- and down-type fermions 
and assure cancellation of anomalies, predicts the existence of three neutral and two charged 
Higgs bosons ho, H, A and H"^. The pursuit of the Higgs bosons predicted in the SM and 
the MSSM is one of the primary goals of the present and next generation of colliders. The 
neutral Higgs bosons productions were intensively studied in the literature [Q. But there are 
relatively less studies on the charged Higgs bosons production. Recently, the calculations of 
the charged Higgs boson pair production at hadron colliders through the gluon-gluon fusion 
have been carried out in Ref. 0, which show that this production mechanism via quark 
and scalar quarks loop would dominate the usual qq annihilation Drell-Yan process if there 
exist sufficiently heavy quarks or scalar quarks, and the production cross section can reach 
a few fb. However, the Next Linear Collider (NLC) operating at a center-of-mass energy of 
500 — 2000GeV^ with the luminosity of the order of lO^^cm^^s^^ can provide an ideal place 
to search for the Higgs boson, especially, it may produce a charged Higgs boson pair with 
much larger observable production rate than at the hadron colliders, because the production 
process can occur at the tree level and isn't suppressed by the Yukawa couplings between the 
light quarks and the Higgs bosons, and the events would be much cleaner than at the hadron 
colliders, meanwhile the parameters of the charged Higgs boson would be easier to extracted. 
In Ref. P], the process e"'"e" -^ H^H~ and the fermion and sfermion one-loop corrections 
to the process at the NLC have been studied, and it pointed out that the corrections are 
typically around —10% in a wide range of the MSSM parameters. Although the tree level 
process e^e~ — >• 77 — > H^H was also already studied almost a fifteen years ago |Q, the 
one-loop radiative corrections to the process are studied just recently [Q. But in Ref. 
the authors considered only the 0{am1/m'^l) Yukawa corrections in the two-Higgs-doublets 
model (2HDM), in which the virtual particles are top and bottom quarks, and the complete 



one-loop radiative corrections to the process hasn't been calculated so far in both of the 
2HDM and the MSSM. 

In this paper, we present the calculations of the squarks loop corrections to the process 
e^e^ ^ ^1^1 ^ H^H~ in the MSSM. Since, in general, the heavy superparticles decouple in 
loop contributions, the contributions from the lighter stop may be important, which would 
lead to observable effect. The complete one- loop supersymmetric corrections to the process 
should include all genuine supersymmetric virtual particles, which will be discussed in details 
elsewhere |^. The structure of this paper is as follows. In Sec. II we give the analytical 
results in terms of the well-known standard notation of the one-loop integrals [^. In Sec. 
Ill we present some numerical examples with discussions. The lengthy expressions of the 
form factors in the amplitude are summarized in Appendix. 

II. CALCULATIONS 

The Feynman diagrams for the process 

7(Pi, Ai) + 7^2, A2) -^ H+{h) + H-ih), (1) 

which include the squarks loop corrections, are shown in Fig. 1-4, where Ai^2 denote the 
helicities of photons. The tree-level diagrams are shown in Fig. |I|. The self-energy, vertex 
and box correction diagrams are depicted in Fig. |^, ^ and ^, respectively. The relevant 
Feynman rules can be found in Ref fl]]. For simplicity, we define the Mandelstam variables 
as 

s = {Pi+P2f = {ki + k2f 
i = (pi - hf = {p2 - k2f 
u = {pi- k2f = {P2 - ki)\ (2) 

and introduce the explicit polarization vectors of the helicities (A1A2) for photons as follows 

e^(pi,Ai = ±l) = --^(0,1,^^,0) 
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e^(p2,A2 = ±l) = -^(0,l,±^,0). (3) 

This choice of polarization vectors for the photons iniphes 

e,.p, = 0, (z,j = l,2) (4) 

and by the momentum conservation, 

ei.ki = -ei.k2, {i = 1,2). (5) 

In our calculations, we used dimensional regularization to control all the ultraviolet 
divergences in the virtual loop corrections and we adopted the on-mass-shell renormalization 
scheme P,^. Note that dimensional reduction is usually more used in the calculations of 
the radiative corrections in the MSSM as it automatically preserves supersymmetry (at least 
to one- loop order), while dimensional regularization does not. But there are only virtual 
scalar particle in the loops in this work, and the calculations of Feynman integrals do not 
involve the n-dimensional Dirac algebra, so the dimensional reduction scheme is same with 
the conventional dimensional regularization one here. Including the squarks loop corrections, 
the renormalized amplitude for the process 77 —>■ H^H~ can be written as 

where Mq is the amplitude at the tree level, (5M'"^'-^, §M^'^'^^'^^ and 6M^°^ represent the 
squarks loop corrections arising from the self-energy, vertex and box diagrams, respectively. 
The corresponding amplitude at the tree level for 77 -^ H^H~ is given by: 

Mo = M* + M^ + Ml (7) 

where Mg, M^, Mq represent the amplitude of the t-channel, u-channel and quadratic cou- 
pling diagrams in Fig. |I|, respectively. Their explicit expressions can be given by 

mH± ^ — t 



rufji — u 
M^ = e\i.e2. (10) 

The amplitude 5M^'^^^ arising from the self-energy corrections is given by 

-Ms ^'":^^:'-^^ ,,z„.,. (u) 

u — mjj± 
where J2^ ^ > 5niH± and 5Zh± represent the charged Higgs self-energy, the charged Higgs 
boson mass renormalization constant and the charged Higgs boson wave function renomal- 
ization constant, respectively, which are given by 

- ^^7 m,-/ 5,(0, m,-/, m,-/) - ^G ^C MO, m,-/, m,-/) 

- i^5 mf^^ Bo{0, mil^, mf^^) -i^e mf/ 5,(0, mj/, m^/) 

+ 6' 5.(A;', m,~/, m,~/) + 6' 5.(A;^ mg/, m,~/)) , (12) 

TT+ TT + 

5ml± = ReJ2 {ml±) (13) 



5Zh± 



\k'^=m? 



+ G' Go{m]j±, mr 2, m^/) + 6' ^^(m^i, m^~/, m,-/)) (14) 



with 



dBo{k^,Mf,M? 



Go{M\ Mf, Mi) = '''''' " W=M^. (15) 



Here we only explicitly write down the third generation scalar quarks contributions, and Nc 
is the number of colors, ^i{i = 1 — 4) and ^i{i = 5 — 8) are the coupling constants of the 



verteces of H^ — t — h and H^ — H — q — q, respectively, and the explicit definitions of 
^i{i = 1,8) are given in Appendix. Bq is the two-point function, definition for which can be 
found in Ref. 0J|. 

The amplitude SM^'^^^'^^ arising from the vertex corrections is given by 

+ bz^A + T^f^^ZzA + ^Z,.) (16) 

2 4 sm dw cos b'ly 

with 



19E|;V1| singH^ E^^(O) ..„s 

^^^ - 2 dk^ '^•^=° ~ cos^vF M| ' ^^^^ 

bZAA = ^p |fc2=0, (isj 

,Z,,^2?U), (19) 

where the expressions of M^^* and Mg" are explicitly presented in the Appendix. From our 
calculations, we find that 

^ZzA = (20) 

for squrks loop diagrams, and the electron charge renormalization constant SZ^. can be can- 
celed by the photon wavefunction renormalization constant ^5Zaa, and only 5Zh± remains 
among the renormalization constants in Eqs. [T^. 

The amplitude 6M^°^ arising from the box corrections is given by 

SM^""" = ei.eaMf"^ + fci.eifci.eaM^^ + 2e2ei.e2<5ZH±, (21) 

where the explicit expressions of Mf°^ and M2°^ are also given in Appendix. Because of the 
same reason as above, we have not written out the SZza, SZaa and 6Ze in counterterms in 
Eqs. 0. 

The corresponding amplitude squared for the process 77 -^ H^H can be written as 



Y, \Mren\ = Y. I^ol + 2ReY{5M'^^f + ,5M^"''*"^ + 5M''°^)MJ, (22) 

where the bar over the summation recalls average over initial photons spins. The cross 
section of process 77 — > H^H~ is 

1 



"max 



^= L TT-^ElMl'di (23) 



^Q^^\pins 



with 

S _ (' 2 "^\ r-^ 

s i3s 

where /3 = Jl — Am?^±/s. The total cross section of e+e~ — >• 77 ^ H^H^ can be obtained 
by folding the a with the photon luminosity 



a{s) = I dz'^-^a{-i-i -^ H^H-at s = zh), (25) 

J2mg±/^ dZ 

where ^/s and v^ is the CMS energy of e+e~ and 77, respectively, and dL^^/dz is the 
photon luminosity, which is defined as 

-^ = 2^y_^ —f,/,ix)f.,/,izyx), (26) 

where f^/e{x) is the photon structure function of the electron beam [l^. For a TeV collider 



with ax/cTy = 25.5, the beamstrahlung photon structure function can be represented as jTl 



'2.25-.,/^) fi^)^^^ for x< 0.84, 
f,/e{x) = \ ^ V 0.166 A ^ ; (27) 

for X > 0.84, 

where x is the relative momentum of the radiated photon and the parent electron. 

If we operate NLC as a mother machine of photon-photon collider in Compton back- 
scattering photon fusion mode, the energy spectrum of the photons is given by |12| 



, T^ (l-x+-r^ ^ + Jf ,^ \ for X < 0.83, Xi = 2(1 + ^2), 

I for X > 0.83. 



III. NUMERICAL EXAMPLES AND CONCLUSIONS 

In the following we present some numerical results for the squarks loop corrections to the 
cross sections of a charged Higgs boson pair production in the process of 77 -^ H^H^ and 
e^e" — >^ 77 — > H^H^, respectively. In our numerical calculations, for the SM parameters, 
we choose m^ = 80.33GeV, m^ = 91.187Ge\^, nit = 176GeV and a = ■^, and other quark 
masses are chosen as zero. Other parameters are determined as follows 

(i) The Higgs boson masses rriho, ^^h-, ^^Ai ^^h± and parameters a, /3 are not constrained 
in the general 2HDM, but in the MSSM, relations |I[] among these parameters are required 
by SUSY, leaving only two parameters free, e.g., mH± and tan/?. Also, in the MSSM the 
charged Higgs boson mass is heavier than the W mass due to the relation 'm'\j± = m^^ + m\. 
In our numerical calculations, we will use mH± as a variable > 100 GeV. 

(ii) In the MSSM the mass eigenstates gi and g2 of the squarks are related to the current 
eigenstates ql and qr by ||l| 






w 



(a ^ 



with R" 



cos^g sm6g ' 



V 



sin 6*5 cos 6*5 



(29) 



For the squarks, the mixing angle 9q and the masses rrig^ 2 can be calculated by diagonizing 
the following mass matrices 






^ M|^ m.MLR^ 



y rrigMRL 



MIr J 



MIl = ml + m^„ + ml cos 2/3(/f - e, sin^ ^^), 



MIr 



WLfj ^ + m^ + ml cos 2/5e„ sin^ 6. 



M, 



LR 



M, 



RL 



q • ■■"z — --^'-g• 

At — fi cot P {(j = i) 
Ah — /itan/? {q = b), 



(30) 



where m^, mj- ^ are soft SUSY breaking mass terms of the left- and right-handed squark, 
respectively. Also, fi is the coefficient of the H1H2 mixing term in the superpotential. At 
and Ab are the coefficient of the dimension-three trilinear soft SUSY-breaking term. J^^, Cq 



8 



are the weak isospin and electric charge of the squark q. From Eqs. ^ and ^, mi^ ^ and % 
can be derived as 

1 



m? = - 

tl,2 2 



ML + ML T J {ML - M|^)2 + 4m|M 



2 
LR 



tan^i = ^^— — ^. (31) 

In our numerical calculations, we always assume that nijj = m^ = rriQ and M^^ > 0. 

The currently most popular supersymmetric models permit to have large splitting in the 
masses of the left-top squarks and right-top squarks, while the masses of all the other (left- 
and right- ) squarks are about the same [|l^. Assuming one has a supersymmetric signal, 
we use the full (> 100 parameters) parameter space freedom of the MSSM and fit the data. 



This approach has been used in studying the CDF e"'"e~77 -|- ^t event [|T4|. It was found 
that to also explain all the low energy data (including a^, Rf, and the branching ratio of 
b -^ S7, etc.), the lightest mass eigenstate (ti) of the top squarks is about to be the order 



of mw', tan/5 is the order of 1; the sign (fi) is negative, and |/i| ~ rriz [0- We shall refer to 
this class of models as MSSM models with scenarios motivated by current data. 

Figure ^ Fig. |^ show the relative corrections 5cr/ao as a function of mH± with vS = 500, 
1000 and 2000 GeV, respectively, assuming m^ = mj^ = ttIq = 1 TeV. The corrections in 
the stops mixing case can become much large when tan/? = 1.5, which can range between 
-20% and 13%, -20% and 25%, -16% and 23% for three different photon-photon CMS 
energy, respectively, but the corrections are all negligibly small in the no-mixing case. In 
order to study the variations of the corrections with respect to the tan j3, we also show the 
curves for tan/5 = 4 and 40, and find that the corrections are much smaller than ones for 
tan/? = 1.5, which are at most a few percent. Here and below, we take the mass of lighter 
stop to be equal to 90 GeV in the mixing case of stops, which is in agreement with current 



experiment lower limit to the squark mass W^. 



In our numerical calculations, we also considered the cases for the small splitting in the 
masses of stops, and found that the corrections are generally small except near threshold, 
which opens the charged Higgs boson decay into two squarks. We don't show the figures 



here for simplicity. 

Fig. ^ gives the total cross sections as a function of charged Higgs boson masses at 
tree-level for ^/s = 500, 1000 and 20006*61^ and two sources of photons modes. For 'mfj± = 
150GeV, the total cross sections can be greater than 100/6 for laser back-scattering photons 
mode whether ^ = 500, lOOOGeF, or 2000Ge\/, and for m^i = 300GeV", the total cross 
sections are a few tens fb and for ^/s = 1000 and 2000 GeV. But for the beamstrahlung 
photons mode, the total cross sections can also reach hundreds of fb for the lighter charged 
Higgs boson masses for ^/s = 1000 and 2000 GeV, and they decrease rapidly with increasing 
the charged Higgs boson masses. 

Fig. H - Fig. present the relative corrections to the total cross sections as a function 
of mfj± with y/s = 500, 1000 and 20006*6^, respectively, assuming m^ = m^, = m^ = 1 
TeV, for two sources of photon modes. In general, those corrections always decrease the 
tree-level total cross sections in the mixing case of stops, which can exceed —10% for a wide 
range of charged Higgs boson mass for tan/? = 1.5. However in other cases, the corrections 
are at most only a few percent reduction. 

In conclusion, we have calculated the squarks one-loop corrections to the cross sections 
of the charged Higgs boson pair production in photon-photon collisions in the MSSM. In 
general, in case of the large splitting in the masses of left-top squark and the right-top 
squark and for tan/3 near 1.5, the corrections can reduce the cross sections by more than 
10% for a wide range of the charged Higgs boson mass, depending on the CMS energy y/s. 
But in other cases, the corrections are at most only a few percent. Those corrections can be 
comparable to the 0{am^/m'^) Yukawa corrections ^. 
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APPENDIX 

In this appendix, we present the exphcit analytic expressions of the form factors only 
including the third generation squarks contributions. 
Ml" is given by 

12 



M*" = E/i 



{tri,i) 



(32) 



i=l 



f{tri,i) 



where the form factor /| ' represents the contributions arising from Feynman diagrams 
with the indices of i in Fig. 3. Here 



n 



(trij) 



for j = 1,2,3,4,5,6,7,8,11,12, 



(33) 



/i 



(tri,9) ie^grriwiT. cos(q; — /?) — cos(2/9) cos(a + P)) 



47r2(s — nilj) 
66efeCoo(0, 0, s, ml^,ml^,ml^ 



+ 
+ 
+ 

+ 
+ 



66 ^^15, mi 



b, - ^6i 



66 -^ et,^i6 -^6o,"^6 



h2 



%2 - "^ti 



ej, ^ et.iiQ -^ i9,mi^ 

igmyj{2 sin(Q; — /?) — cos(2/3) sin(Q; + /3)) 

47r2(m^o - s) 
69efeCoo(0, 0, s, ml^,ml^,ml^] 



69 -^Gs,"^, 



fe ^ "^61 



66 ^ et, ^9 -^ 63, "^6 



fe 



66 ^64,^9 ^62, "^6 



62 ^ "^ti 



"^t2 



"^t2 



(34) 



/i 



trj.io ie^gmu,{2 cos[a — /?) — cos(2/5) cos(« + /?)) 



167r^(m|^ — s) 



(666550(5, 



<'<^ 



+ 6565^0(5, m?^,m?J + ^ioejSo(s,m|,m|) + ^Qe^Bo{s,ml,ml 



igmw{2 sin(Q; — /9) — cos(2/9) sin(a + /?)) 



167r2(s - m^o) 



(ei96^5o( 



+ 6866^0(5, m^^,m^ J + ^i3et5o(s,mj^,m^-J + ^126^^0 (s,m^-^,mj J j , 



(35) 
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where et, Cb are the electric charges of the top and bottom quarks, respectively. Here and 
below, Bjn, Cm, Cmn, -Dm, Djnn are the two-, three-, four-point Feynman integrals, definitions 
of which can be found in Ref. |0,§]. C,j, (j = 9, 10) is the coupling constants of the vertex 
H — t — t, which are given by 



?9 Ul 



en e 
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{-zg)R\A^i,){R') 



i\T 



with 



^fTfi 



I m^((l/2)-(2/3)sin^6>,„)cos(a+/3) , m^ sin(a) 



Hit 



mt 



mm sin(/3) 2™,^, sin/3 



{At sin a — ^ cos a) 



\ 



V 



mt 



2miu sin (9 



^-g {At sin a — /i cos a) 



mz(2/3) sin^ 9w cos(a+/3) _|_ r>i| sin a 

cos 9u, m,„ sin /3 / 



^j, (j = 12, 13) is the coupling constants of the vertex Hq —t — t, which are given by 



^12 ^4 



V 



u e 



13 



with 



A 



hoii 



( m.((l/2)-(2/3)sin^e.)sin(a+/3) _ rr^^^ _r2^(_^^ cOS « - /i siu a^ ^ 

mz (2/3) sin^ 6>m sin(a+/3) m^ cos a 



\ 



2mu, sin/3 ^ 



-y4( COS a — yU sin a) 



^/ 



0' (j 



15, 16) is the coupling constants of the vertex H — b — b, which are given by 

H9)R\A^-J{Rr 



^15 ^17 



in i 



16 



with 



/ -mz((l/2)-(l/3) sin^ d^) cos(a+/3) , mgcos(a) 



A 



Hbb 



+ 



m,, 



V 



cos6,u muicos(p) 2mwCosf3 

mt ( A „^„ ^, ,, „• ^\ -mz „;„2 



(AfjCosa — /isina) 



\ 



((Afocosa — /isina) 



sin Qyj cos(q; + /?) 



mf cos a 



2mu, cos/3 ^^ "^^ ^^'^ " ^ uiij. i^y 3cos9ti, '^^^'^ v^^, v^vj^jy^ix, i ^y i ra^co&fi I 

C,j, {j = 18, 19) is the coupling constants of the vertex Hq — b — b, which are given by 



U8 ^20 



6o e 



19 



{^g)R\A.^i){R 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



(42) 
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with 



A 



' -m^((l/2)-(l/3)sin^6>^)sin(a+/3) , mj sin(a) 
— ^Q ~r ^^ ( a\ 



hobb 



rrit 



V 



COS 6. 

-^ — j {Ah COS a + fi sin a 



niw cos(/3) 2m,^ cos /3 

2 



(Afe COS a + /i sin a) 



2mw cos /3 



3cos( 



sin 9^ sm{a + (3) + 



■ (43) 



^/ 



M2 * is given by 



12 






(44) 



p(tri,j) 



where the form factor /2 ' represents the contributions arising from Feynman diagrams 



with the indices of i in Fig. 3. Here 



/i*"'^') = 0forj = l,2,3,4,9,10,ll,12, 



(45) 



f: 



(tri,5) 



47r2(m2^±-t) 



^|e6C2(m^±, 0, i, m\,ml^,ml^) - iletC2{rn]j±,Q, i, m\,m\,ml 



b2' 



+ 



^3 -* 6, "^ti -^ rrii^ + 6 ^ ^4, ^ii ^ "^t2, "^; 



^b, -^ ^h 



+ 



6 ^^i,m. 



62 ^ ^61 



, (46) 



f: 



(trifi) 



^3et(Co(0,m^±,t,m^-^,m^-^,m^J + (Ci + C2)(m^±, 0, t,m^^,m^-^,m^-^ 



A-K'^{m]j± -t) 
- ^|e6(Co(0, m^±, t, ml^,ml^,m\) + (C^ + (:72)("^h±5 0, t, m\,ml^,ml^) 



+ 



6^6, "^ti ^ "^t2 + 6 ^ ^4, rrii^ -^ rni^,mh 



b2 - ^61 



6 ^6,"^L 



62 - ^61 



, (47) 



/i*"''^ = /i*"''^(^-t,A;i^A;2) 



(48) 



ft''^ = ft''\u^i,h^h) 



(49) 



where ^j, (j = 1, 2, 3, 4) is the couphng constants of the vertex H^ — t — b, which are given 
by 






i?* 



«^ 



V2^. 



y ^4 6 y 

M^*"^ is given by 



m„ 



/ 



— m^ sin(2/9) + mfcotp + mltanP mb{fi + AbtanP) 
mt{fi + Atcotp) mhmtitanP + cotp) 



\ 



{Ry. (50) 
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M'r = E /i 



(boXyi) 



(51) 



i=l 



f-{box,i) 



where the form factor f{ ' represents the contributions arising from Feynman diagrams 
with the indices of i in Fig. 4. Here 



n 



{box,l) 



47r2 



V2 2 T^ / 2 n n 2 2 ^ 2 2 2 2 



+ ^oo("^H±, 0, 0, m^±, t, s, m^^,m-^^,m-^^,m-^^ 



+ 



6 ^6,"^fei 



m 



&2 



6 ^ ^i,m 



t2 



mr 



+ 



"^t2 -^^h^^h ^^h 



(52) 



/i 



(6ox,2) _ r{box,l) 



fl ' '{t ^u,ki^ k2), 



(53) 



/i 



(box,3) _ e e&gf 



47r2 



([e|^oo(0, 



2 r, 2 -X2 2 2 2^ 

mH± , 0, m^± , u, t, m^^ , m^^ , m^^ , m^^ ^ 



+ Doo(0,m^±,0,m^±,M,t,m^^,m^^,m^-^,m^-J 



+ 



Cs -^6,"^ti 



^2 



^3 ^^i,mi 



b2 



^fel 



mi 



^t2'"^b2 



"^fel 



, (54) 



fl 



{box,4) 



ij^e^elCQo (0, 0, s, m\ , m| , mf^ ] 

27^ 



^h ^"^t2'^5 ^^6 



+ 



^h -^ ^h > ^5 ^ 6 + "^ii ^ "^fe ,6^6 



(55) 



A^ 



(feox,5) 



87r2 



yz,^ /2 2-2 2 2\i yZ/^ I 2 2-2 2 2\ 

^4C;o(m^±, mj;^±, s, m^-^, m^^^,m^J + ^iC;o(m^±, m^±, s, m^^,m-^^,mfj 



I ^^2/^ /2 2-2 2 2\i /•2/^ / 2 2-2 2 2 ' 

+ ^2 Co {mH± , m^± , s, m^-^ , m^^ , m^-^ ) + C3 Q (m^± , m^± , s, m^^ , m^^ , m^^ ^ 



+ 



87r2 



^h ^^h^^h ^^b2 



(56) 



/I 



{box,6) 



—le 



i-K 



2 y^selBois, ml^,ml^) + ^76^^0(5, ml^,ml^) + UelB^i^s, m\,m\] 



+ ^5etBois,m^^,m^^ 



(57) 



where ^j, (j = 5,6,7,8) is the couphng constants of the vertexes H~^ — H — t — t and 
H^ — H — b — b, respectively, which are given by 
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with Xi, {i = 1, 2, 3, 4) is irrelevant to our calculations. 
The form factor Mf""^ is given by 



M. 



box 
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(60) 



where the form factor /2 °^' represents the contributions arising from Feynman diagrams 
with the indices of i in Fig. 4. Here 
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fr'" = for J = 4, 5, 6. (64) 



[1] H.E. Haber and G.L. Kane, Phys. Rep. 117, 75(1985); J.F. Gunion and H.E. Haber, Nucl. 
Phys. B272, 1(1986); For a review see, e.g. J.F. Gunion, Higgs Hunter's Guide (Addison- 
Wesley, Reading, MA, 1990). 

[2] Y. Jiang et.al., J.Phys. G 23, 385 (1997); Erratum-ibid, J. Phys. G 23, 1151 (1997); Y. Jinag 
et.al.. Report No hep-ph/9708421; S.S.D. Willenbrock, Phys. Rev. D 35, 173 (1987); A. Krause, 
T. Plehn, M. Spira and P.M. Zerwas, Report No. hep-ph/9707430. 

[3] A. Arhrib, M. C. Peyranere, and G. Moultaka, Phys. Lett. B341,313 (1995). 

[4] D. Bowser-Chao, K. Cheung and S. Thomas, Phys. Lett. B 315, 399 (1983); M. Drees, R.M. 
Godbole, M. Nowakowski and S. D. Rindani, Report No. MAD/PH/824. 

[5] W.G. Ma, C.S. Li and L. Han, Phys. Rev. D53, 1304 (1996); Erratum-ibid, Phys. Rev. D 56, 
4420 (1997). 

[6] C.S. Li, S.H. Zhu, C.S. Gao, in preparation. 

[7] G. Passarino and M. Veltman, Nucl. Phys. B 33, 151 (1979); R. Mertig, M. Bohm and A. 
Denner, Comp. Phys. Comm. 64, 345 (1991). 

[8] A. Denner, Fortschr. Phys. 41, 307 (1993). 

[9] S. Sirlin, Phys. Rev. D 22, 971 (1980); W. J. Marciano and A. Sirhn, ibid. 22, 2695 (1980); 31, 
213(E) (1985); A. Sirhn and W.J. Marciano, Nucl. Phys. B 189, 442 (1981); K.I. Aoki et.al.. 
Prog. Theor. Phys. Suppl. 73, 1 (1982). 

[10] F.Halzen, C.S.Kim and M.L.Stong, Phys. Lett. B274, 489(1992). 

[11] R. Blankenbecler and S.D.Drell, Phys. Rev. Lett. 61, 2324 (1988); D.V.Schroeder, SLAC- 



16 



Report-371 (1990). 

[12] M.Bohm, W.Hollik and H. Spiesberger, Fortschr. Phys. C27, 523 (1985). 

[13] For example, see G.L. Kane, C Kolda, L. Roszkowski, J.D. Wells, Phys. Rev. D49 (1994)6173, 
and references therein. 

[14] S.Ambrosanio, G.L. Kane, G.D. Kribs, S.P. Martin and S. Mrenna, Phys. Rev. Lett. 
76(1996)3498; S.Dimopoulos, M. Dine, S. Raby and S. Thomas, Phys.Rev.Lett. 76(1996)3494. 

[15] R. Barate et.al., Phys.Lett. B413 (1997) 431. 

FIGURES CAPTION 

Figure [^: The tree level diagrams for sub-process 77 -^ H^H^ . In the Feynman di- 
agrams of this figure and below, the wavy, the dashed and the dotted lines stand for the 
photon, the charged Higgs boson and the squark, respectively. 

Figure 0: The self-energy correction diagrams of sub-process 77 -^ H^H . 

Figure |^: The vertex correction diagrams of sub-process 77 — > H^H^ . 

Figure |^: The box correction diagrams of sub-process 77 -^ H^H . 

Figure |^: Relative corrections as a function of the charged Higgs boson mass with ys = 
SOOGeV, where no-mixing represents that there does not exist mixing between the left- and 
right-handed stop, and mixing represents there is mixing between the stops. We choose 
rrifj = nijj = ttIq = 1 TeV and fi = —90 GeV. In the case of mixing of stops, we choose the 
value of parameter At which makes the masse of the lighter stop equal to 90 GeV. 

Figure ||: Same with Figure | but Vs = lOOOGeV^. 

Figure 0: Same with Figure | but v^ = 2000GeV^. 

Figure ^. The tree level cross sections for the process e^e^ ^ 77 ^ H^H^ as func- 
tions of the charged Higgs boson masses, where the e^e~ CMS energies are 2000, 1000 
and SOOGeV^, respectively. The letters (L) and (B) on the curves represent the laser-back 
scattering and beamstrahlung photons modes. 



Figure |^: Relative corrections as a function of the charged Higgs boson mass with ^/s = 
500GeV, where no-mixing represents that there does not exist mixing between the left- and 
right-handed stop, and mixing represents there is mixing between the stops. We choose 



ruf 



m 



f) = uIq = 1 TeV and fi = —90 GeV. In the case of mixing of stops, we choose the 



value of parameter At which makes the masse of the lighter stop equal to 90GeV. The letters 
(L) and (B) on the curves represent the laser-back scattering and beamstrahlung photons 
modes. 

Figure [l^: Same with Figure || but the e~^e~ CMS energy y/s = lOOOGeV^. 

Figure [Tl|: Same with Figure || but the e"'"e~ CMS energy ^i = 20006*61/^. 
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